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SUMMARY 

The nonequilibrium chemical reaction of dis- 
sociation and recombination is studied theoretically 
for air in the viscous shock layer at the stagnation 
region of axisymmetric bodies. The flight regime 
considered is for speeds near satellite speed and for 
altitudes between 200,000 and 300,000 feet. 

As altitude is increased from 200,000 to 300,000 
feel , the flow in the shock layer varies from one which 
supports a well-defined boundary layer to one 
characterized by a thick viscous layer. The chemical 
reaction rate , at the same time , varies from near 
equilibrium to completely frozen. 

The nonequilibrium chemical reaction is found to 
produce a significant reduction in the heat transfer 
to none at aly tic walls at certain altitudes. 

It is shown that the controlling chemical reaction 
for heat transfer at higher altitudes is the dissociative 
shock-layer reaction, whereas at lower altitudes it is 
the bound ary -layer recombination process. 

The effects of nose radius and wall temperature on 
the heat transfer are studied. 

The flight condition at which the air at the bound- 
ary-layer edge reaches a state of chemical equilibrium 
is found 

Also , the simultaneous effect of air rarefaction 
and nonequilibrium chemical reaction on the shock 
layer thickness is analyzed. 

INTRODUCTION 

The flow of chemically reacting gases has been a 
subject of considerable interest in recent years. 
The subject became particularly important in 
external aerodynamics in connection with the 
calculation of skin friction and heat transfer for 
hypersonic vehicles. The classical boundary- 
layer theory usually can be used to analyze skin 
friction and heat transfer for flights of blunt bodies 


at altitudes below about 200,000 feet (refs. 1, 2, 
and 3). For such a flight regime, chemical re- 
action processes in the in viscid region and in the 
boundary layer can be studied independently. 
Studies of the nonequilibrium dissociation and 
recombination have been carried out by various 
investigators for the in viscid region (refs. 4 through 
6) and for the boundary layer (refs. 7 through 9). 

With the current interest in manned space flight 
and entry into the earth's atmosphere, study of 
nonequilibrium aerodynamics must be extended 
beyond 200,000 feet altitude. This is true because 
the maximum convective beating for many blunt- 
nosed vehicles returning to the earth from para- 
bolic orbits may occur at altitudes above 200,000 
feet (ref. 10). In addition to the problem of con- 
vective heat transfer, an understanding of the 
chemical reaction in the entire shock layer near the 
stagnation region is very important in the analysis 
of radiative heat transfer, communication, and 
re-entry trails (refs. 11 and J2). The study of the 
nonequilibrium viscous shock layer becomes in- 
creasingly complicated as altitudes are increased 
above about 200,000 feet because the in viscid and 
the viscous flow regimes can no longer he anahzed 
separately. The complication is due to the inter- 
action of vorticity and chemical reaction in the 
zone between the shock layer and the boundary 
layer. 

The viscous shock layer in the rarefied atmos- 
phere between the altitudes of 200,000 feet and 
300,000 feet was studied in references 1, 2, and 3 
with the air assumed to be chemically inert. The 
analysis showed that, at the stagnation region of 
spheres, the skin friction and heat transfer are 
much higher than those calculated by classical 
boundary-layer theory in the range of 10 3 ^ Re ^ 10 5 . 
These Reynolds numbers roughly correspond to 
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the flight regime of present interest for blunt bodies 
with nose radii between 1 and 5 feet. 

In the present paper, the chemical reaction of 
nonequilibrium dissociation and recombination in 
the viscous shock layer will be studied, and the 
effect of the chemical reaction on stagnation point 
heat transfer will be investigated for noncatalytic 
walls. The flight regime of the blunt bodies con- 
sidered in the present study is for altitudes between 
200,000 feet and 300,000 feet and for flight speeds 
near satellite speed. 

As will be seen in the text, the actual numerical 
solutions of the governing equations require con- 
siderable care because of the complications im- 
posed by the nonequilibrium effect. In order not 
to further complicate the numerical phase of the 
analysis, certain simplifying assumptions will be 
made on the fluid properties. 

SYMBOLS 

a speed of sound 

B constant defined by equation (A8) 

C atom mass fraction 

( \ mass fraction of ith species 

(\ complementary solution of equation (Al) 
and solution of equation (A4) 

(\ particular solution of equation (Al) and 
solution of equation (A5) 

K n energy of the nth electronic state 
f dimensionless stream function defined by 
equations (41) and (42) 
g n degeneracy of the rttli electronic state 



h static enthalpy of air including chemical 
energy (~A* in the shock layer) 
hi static enthalpy of ith species 

/ti° heat of formation of ith species 
A A° heat of dissociation 

h t total enthalpy defined by equation (15) 

h Plank constant 

I molecular moment of inertia 

j 0 for two-dimensional bodies; 1 for axisym- 
i lie trie bodies 

K curvature of the body surface 
K p equilibrium constant based on partial 
pressures 

k Boltzmann constant 

k T specific recombination coefficient 

L nose radius 

l characteristic length 


M Mach number 

Mi molecular weight of ith species 
mass of ith particle 
N total number of particles 

n 1,2, 3,. . . 

p pressure 

Pr Prandtl number 

Q t total partition function of ith species 
Q i e electronic partition function of ith species 

Q i r rotational partition function of ith species 

Q it translational partition function of ith 

species 

Q iv vibrational partition function of ith species 
q heat-transfer rate to wall per unit area 
It universal gas constant 

Be free-stream Reynolds number, Pa>Uao ^ J 

T absolute temperature 

u x component of velocity 

hi x component of velocity gradient 

u m free-stream velocity 

v y component of velocity 

W production rate of atoms by chemical re- 

action per unit volume 
A r ( mole concentration of ith species 

x direction parallel to wall (see fig. 1) 

y direction normal to wall (see fig. 1) 

ft, streainwise velocity gradient at the 

boundary-layer edge 
T constant used in the appendix 

A shock-layer thickness 

e density ratio across shock, — - 

Ps 

r] similarity variable defined by equation (40) 

0 D characteristic temperature of dissociation 

6 V characteristic vibrational temperature 

A mean free path 

H dynamic, coefficient of viscosity 

/ dilat at ional coefficient of viscosity 

g" longitudinal coefficient of viscosity 

p density 

SUBSCRIPTS 

AK adiabatic equilibrium, the value which 

would be obtained if the air were allowed 
to reach equilibrium adiahaticallv with 
K 

BL boundary layer 

K local equilibrium 

i ith species 

n nth state or nth iteration 

immediately behind shock 
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SL shock layer 

w wall 

! atoms 

2 molecules 

3 third body which is either atom or molecule 

co free stream 

SUPERSCRIPT 

' total differentiation with respect to rj 

SHOCK-LAYER CONSERVATION EQUATIONS 

In the present analysis, the simplified forms of 
the Navier -Stokes equation and the corresponding 
energy and diffusion equations will he used in con- 
junction with the no-slip wall conditions and the 
conventional shock relations. We shall first try 
to justify the use of the equations and the wall and 
shock conditions for the flight regime of altitudes 
between 200,000 feet and .‘100,000 leet and for 
Mach numbers of near 30. We shall then develop 
the appropriate conservation equations and bound- 
ary conditions. 

In 1940, Tsien (ref. 13) first proposed a classifi- 
cation of fluid mechanics based on (he Knudscn 
number. Knudsen number expresses the degree 
of rarefaction of the fluid and is defined as the ratio 
of the mean free path to a characteristic length, 
\/L Of the various regimes into which Tsien and 
others catagorized gasdynamics, in general, only 
two are reasonably well defined. They are the 
regime of ordinary continuum gasdynamics where- 
in the density is sufficiently high that inter- 
molecular collisions dominate over collisions with 
the boundaries, and the regime of free molecule 
flow where the gas is sufficiently rarefied that 
collisions with the boundaries dominate over col- 
isions between the molecules. In terms of the 
Knudsen number, the two regions arc characterized 
by the inequalities A//<C<0 and 1, 

respectively. 

Since 1946, considerable theoretical work has 
been done on the subject of free molecule flow 
(see, e.g., ref. 14). The problem of free molecule 
flow was first attacked because the classical kinetic 
theory could be applied to this limiting case more 
readily than to the intermediate regimes. 

In recent years, the intermediate regime has 
been studied with new interest and, as a conse- 
quence, some promising theories have been formu- 
lated and a few interesting solutions have been 
obtained. Some of this work appears in references 
1, 2, 3, and 15. The recent success in studying the 


intermediate range of rarefied gasdynamics was 
made possible by the continuum approach which 
utilizes the Xavier-Stokes equation. Probstein 
and Kemp in reference 1 subdivided the entire 
flow regime, between t lie boundary-layer regime 
and the free molecule regime, into seven sub- 
regimes: boundary -layer regime, vorticity inter- 

action regime, viscous-layer regime, incipient 
merged layer regime, fully merged layer regime, 
first-order collision theory regime, and free mole- 
cule flow regime. Probstein and Kemp suggested 
that the Navier-Stokes equation be used through 
the fully merged layer regime for engineering 
purposes. Moreover, through the viscous-layer 
regime, the shock wave near the stagnation region 
may be considered as a finite discontinuity and the 
conventional shock relationships may be used to 
predict the flow conditions immediately behind 
the shock wave. The justifications for these sug- 
gestions are found in references 1 and 15. We shall 
here briefly recapitulate the portions of the argu- 
ments which are directly pertinent to the present 
work. 

We shall first show that the flight regime of 
present interest lies in the continuum and the 
viscous-layer regimes. In the present paper, we 
are interested in the shock layer at the stagnation 
region of axisymmetric bodies only. The mean 
free path and the characteristic length ol concern 
are, then, those in the shock layer itself. Tn order 
that the continuum theory be valid in the present 
study, we must have (W^)«l. The mean free 
path \ (SL is longest immediately behind the shock 
wave where the effect of the cooled wall is not 
felt. We shall, therefore, justify the continuum 
theory on the basis of the mean free path immedi- 
ately behind the shock to be on the safe side. 
The continuum theory is then valid when 

a) 

With the aid of kinetic theory, we have 



( 2 ) 

(3) 


With the aid of equations (2) and (3), the con- 
tinuum criterion (1) becomes 
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i&«‘ < 4 > 

The conventional shock relations can be used, in 
general, if the conductive process expressed by 
f*s{bu/dy) s is much smaller than the convective 
process given by p s v s v t at the shock wave. With 
the aid of the assumptions (du/dy )^^/ A and 
ttcoi it can be shown that the magnitude 
of the conductive process becomes much smaller 
than that of the convective process when 

Mm ^ ^ r { r\ 

Re «' v t 

Now consider a body with a nose radius of 1 foot 
traveling at the speed of = At the alti- 

tudes of 250,000 and 300,000 feet, the free-streain 
Reynolds numbers are in the order of 10 4 and 10 3 , 
respectively. The ratio MJIte is, therefore, of 
the orders of 0.003 and 0.03, respectively. The 
density ratio e is the order of 0.1 at these flight 
conditions. It is seen, therefore, that the con- 
tinuum condition (4) and the viscous layer con- 
dition (5) are satisfied at the stagnation point of 
axis vm metric bodies for flight altitudes below 
about 300,000 feet. 

Next we shall show that the satisfaction of the 
continuum criterion (4) automatically validates 
the no-slip wall conditions for a cooled wall. The 
velocity jump at the wall can be expressed as: 



The effect of slip will be negligible in the solutions 
of the present problem if 

(7) 

At the rarefied condition where slip may begin to 
occur, the shear stress p(du/dy) will be of the 
same order of magnitude all along the stagnation 
stream line from the body to the shock. Inequal- 
ity (7) becomes, with the aid of the relationships 
(2), (3), and (6) 



For the present case of T w «T Si yiT s /T a >l and, 
therefore, condition (8) is automatically met when 
the continuum criterion (4) is satisfied. It is seen, 


theref >re, that slip is not an important phenom- 
enon in the stagnation region in affecting skin 
friction and heat transfer. 

Now it is seen that, for the flight regime of 
present interest, we may use the Navier-Stokes 
equation and the corresponding forms of the 
energy and diffusion equations to study the vis- 
cous shock layer at the stagnation region of 
axisymmetric bodies. It is also seen that the no- 
slip wall conditions and the conventional shock 
relationships may be applied to the present case. 

In tne remainder of this section we shall develop 
a set of shock-layer conservation equations. 
Though the development given is brief, the casual 
reader may skip over to the final form of the 
conservation equations, (18) through (28). Aside 
from the nonequilibrium chemical effect, the equa- 
tions used in the subsequent portion are similar to 
those developed in reference 15. 

For a reasonably thin shock layer over two- 
dimensional and axisymmetric bodies, the Navier- 
Stokes, the energy, and the diffusion equations 
can be written in the following simplified form 
for a boundary-hiyer-type coordinate system. 
(See fig. 1 for the coordinate system.) 


Con tii uit.y: 



Figure 1 . — Physical model considered. 
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x component of momentum 

bu . bu ,, 1 dp Id/ bu\ ... n . 


y component of momentum 

u Ku 2 -\~- ~==- 

ox Oy poyp Oy 


«) 


0 2 u 
p bxdy 


Energy 


u 


Diff usion 


ae i d 

U dx + V dy pby\Prby 


(>- 5 >r . ') + 

\Pr by / 


dr\ , tj; 

P 


where 




h, = h+ 2 

=2/’ ( A,+ W ^- 


(ID 


(12) 


(i») 

(14) 

(15) 


Lewis number is assumed to be unity in the 
energy and diffusion equations. In the deriva- 
tion of the diffusion equation (13), the air is 
assumed to be a molecule-atom binary mixture 
and tlie effect of thermal diffusion is neglected. 
The justification for neglecting the thermal 
diffusion effect is as follows. First of all, the 
effect of thermal diffusion is negligible compared 
to that of concentration diffusion at temperatures 
below about 10,000° K (refs. 7 and 15). The 
effect of thermal diffusion is, therefore, not 
important near the cooled wall. Near the shock 
wave, on the other hand, the temperature may 
be considerably higher than 10,000° Iv. The 
convective process, however, predominates over 
the conductive process near the shock wave 
since condition (5) is met in the present flight 
regime. Thermal diffusion is, therefore, neglected 
in the present study. 

Now, for the purpose of simplicity, we shall 
make the assumption that the shock laj^er is 
very thin, so that 



(16) 


The above assumption makes the results of the 
subsequent analysis quite restrictive. The as- 
sumption, however, should not cause any undue 
errors in the present analysis which is concerned 
with extremely large Mach numbers only. The 
following order-of-magnitude arguments are now 
used to simplify the equations further. 


A 

L 


P <x 

p 


v p 
— ~ 




P h t t 

'V/ I 

P^Uj uj 


(17) 


Guided by experience with stagnation-point 
boundary layer work, it is assumed that all the 
fluid properties and the dependent variables, 
except the x component of the velocity and 
pressure, are functions of y only. The velocity 
component u is expressed as u~Ui(y)x y and 
also, the variation of pressure in the x direction 
is expressed by Newtonian theory. The con- 
servation equations (9) through (13) are now 
simplified for a sphere as: 


Continuity 



2pM,+, |/ -0 

(IS) 

Momentum 

„ , (lu,\ 2p 1 (l / diii\ 

(19) 

u 

1 <h p// p<ly\ <ly) 

Energy 

dh l d / p dh\ 

(20) 1 


L dy p dy \Pr dy) 

Diffusion 

d( 7_1 d / p dC\ W 
L dy~ pdy\Prdy / + p 

(21) 

The equation of state is 



,,^p{\+C)~T 

(22) 


The boundary conditions are as follows: At the 
wall (y = 0), 

Mi— r=0 (23) 

h—h w (24) 

-j — — 0 (non catalytic wall) (25) 


1 In the shock layer, (u2+r 2 )/2<<A; therefore we let ft ( =A. 
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at I lie shock wave (//=?/,), 


L 


->v 


PqU co 

p* 


ll 



(26) 


(27) 


C- 0 


(2S) 


As shown by the boundary condition (25), the wall 
is considered to he nonc-atalytic. Also, the chemi- 
cal reaction is considered to be frozen through the 
shock wave and, therefore, the boundary condition 
(28) is valid. 

The conservation equations and the boundary 
conditions, equations (18) through (28), will be- 
come complete and ready for solution as soon as 
t lie expression for the chemical reaction rate term 
W p becomes known. In the next section, the 
gas model and the chemical reaction kinetics will 
be discussed, A simplified expression for the 
enthalpy, A, and the expression for the chemical 
reaction rate IF/p will be. derived. 

GAS MODEL AND CHEMICAL REACTION 

in the present section, the gas model and the 
chemical reaction will he discussed for the air in 
the shock layer. As a consequence, the simpli- 
fied (expression for the enthalpy and an expression 
for the chemical reaction rate will be derived. 

For simplicity, we shall consider that the air in 
the shock layer is a binary mixture of air atoms 
and air molecules. We shall assume that the 
rotational and the vibrational degrees of freedom 
are in equilibrium with the translational degree of 
freedom at all times. All the electronic degrees 
of freedom, on the other hand, are assumed to be 
unexcited, and it is also assumed that ionization 
is negligible. The lioneqiiilibrium phenomena 
considered are, therefore, those of dissociation 
and recombination only. 

The rotational degree of freedom probably is in 
equilibrium with the translational degree of free- 
dom. The vibrational temperature, on the other 
hand, may lag behind the translational tempera- 
ture because of the rather long vibrational relaxa- 
tion time near the shock wave. The vibrational 
relaxation, under such circumstances, may be 
coupled with the dissociative relaxation, and the 
two relaxation phenomena may occur simul- 


taneously near the shock wave. The addition of 
the simultaneous relaxation phenomenon increases 
the complexity of the solution considerably; 
therefore, we shall leave this phenomenon to a 
later study. The effect of finite vibrational re- 
laxation, however, should not alter the present 
results greatly since the thermal effect associated 
with the vibrational mode of freedom is much less 
than that associated with dissociation and recom- 
bination. The electronic states probably will be 
excite 1 to some degree at the elevated temperature 
near the shock wave. The excitation of the elec- 
tronic states affects the present problem basically 
through its effect on the partition functions. A 
study of the partition functions of oxygen and 
nitrogen atoms and molecules showed that the 
contribution of the fully excited electronic states 
to the partition functions is very small compared 
to that of the other degrees of freedom (see ref. 
lf>). The assumption of unexcited electronic 
states may he justified here for engineering pur- 
poses. Finally, there could be some ionization 
of the hot air molecules near the shock wave. 
Most of the ionization, however, will probably 
take place following the dissociation. The tem- 
perature will decrease rapidly with dissociation 
and this will lend to decrease the ionization. At 
present, nothing can he said with any certainty 
on the- amount of ionization which may be present 
at th< high temperatures mainly because of the 
lack o ' information on high temperature ionization 
reaction kinetics. Moreover, the presence of a 
large amount of ionization will invalidate the 
prose i t conservation equations, Jn the present 
study we shall neglect ionization completely 
and, i t the same time, limit the analysis to the 
flight speeds up to about 30,000 ft /sec. 

Now with the use of the above mentioned sim- 
plifications, we shall develop the expressions for 
the ertlialpy and the chemical reaction rate. 

In order to obtain the properties of the air, we 
must irst obtain the total partition functions for 
the atoms and the molecules. Making the usual 
assumption that no coupling exists between the 
different modes of energy, we have, 

Qi = Qi.tQi.rQiMi.e (29) 

The general expressions of the partial partition 
functions can be obtained from the standard 
literature on statistical mechanics (e.g., refs. 16 
and 17). They are: 
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Ql,r=0 

0 4ir‘ 2 IkT 

H- r a 2 

h k 

Ql.r = » 

,z= \—e~ >t,r 


(30) 


The enthalpy of the fluid is obtained by the 
equation 


h=2C i h l ='20 1 



(31) 


For the present atom-molecule binary mixture 
with unexcited electronic states, equation (II 1) 
becomes with the aid of equations (30): 



+5C-; <1-0 


(rfPrftaf-*'* 


(32) 


The term (1 — (> r )[(d r j7 1 )/(e 6v ' T — ])] in equation (32) 
repres( k nts the contribution of the vibrational 
degree of freedom of the molecules to the enthalpy. 
The value of the term varies from 0 to (1— ( r ) as 
the temperature is varied from absolute zero to 
infinity. The contribution of the vibrational de- 
gree of freedom to enthalpy is maximum when 
('— -0. The maximum contribution, according to 
equation (32), is about 20 percent of the total 
enthalpy, provided no abnormal vibrational ex- 
citations exist. The value of the term, (OJT)/ 
(e 0t?/T — 1 ) , becomes practically unity at a tempera- 
ture of about 4000° K. Lighthill, in reference 18, 
set (0jT)/(e*' fT -l) = lf2 because t ho maximum 
temperature considered was not much above 
4000° K. In the flight regime of present interest, 
however, the maximum temperature is much 
higher than 4000° K, and most of the chemical 
reactions take place at temperatures above 4000° 
lv. W o shall, therefore, set (9 v /T)/(e 6l ' /T —\) — \ 
in the present analysis. Equation (32) now 
becomes 

/<=(!+T ; )^n CaA° (33) 


The molecular weight and the heat of dissociation 
of air used in the subsequent numerical work are 
the weighed average values of () 2 and No. 

The type of chemical reaction considered in the 
present paper is: 

Ah+Ah+A^X+A* (34) 


For the chemical reaction expressed by equation 
(34), the applications of the law of mass action 
and the equation of state yield the following 
equation for the reaction rate: 




In order for equation (35) to be quantitatively 
useful, the equilibrium constant K p and the specific 
recombination coefficient. k r must be known. 

The equilibrium constant can be quite accurate- 
ly calculated by the use of the partition functions. 
At chemical equilibrium, the concentrations of 
the species are given by the following equation 
(see refs. 16 and 17): 

(number of atoms)- Qr - 9 /T 

number of molecules (), C 


The equilibrium constant for the 
by (‘(illation (34) is then derived 
(36) as 


K t 


P_ Qr 


f>n t 


reaction shown 
from equation 


(37) 


The specific expanded expressions for the partial 
partition functions, equations (30), are given for 
the atoms aud molecules of oxygen and nitrogen 
in reference 10. The expressions given in reference 
10 are for temperatures up to only 15,000° K. 
With the present assumption of unexcited elec- 
tronic states, however, these expressions can he 
used in the present study even for temperatures 
much above 15,000° K. The specific properties, 
such as d Vy 0 D , etc., used in the subsequent numeri- 
cal work are the weighed average values of oxygen 
and nitrogen. 

The exact values of the specific recombination 
coefficient k r are not so readily obtainable by a 
theoretical treatment alone. All the reliable ex- 
perimental data available so far are limited to 
temperatures up to about 7000° K. Experimental 
results of specific recombination coefficients are 
reported in references HI and 20 for oxygen and 
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nitrogen for temperatures up to about 7000° K. 
The results showed that the recombination coeffi- 
cients varied with the third body which acts as a 
catalyst in the chemical reaction. A study of the 
experimental results of the references and those 
obtained by other investigators showed that the 
specific recombination coefficient for the air in t lie 
shock layer, a mixture of atoms and molecules of 
oxygen and nitrogen, could be represented reason- 
ably well by the equation 


/ 1 \ L ■"> ( 

2X-,=1 ..j(*X10 20 ( 7T , ) — , — 

\1 / mole 2 sec 


(3«s) 


The same recombination coefficient given by the 
above equation was also used in the boundary- 
layer work of reference 9. 

The reaction rate. U p, given by equation (35) 
now becomes, with the aid of equations (37) and 
(38), 


W 

— - ~2(1.5GX 10 20 ) 

p 



1 (jl91 

4 p \x Q: 



where / and }> are in degrees Kelvin and atmos- 
pheres, respectively. 


SOLUTION OF CONSERVATION EQUATIONS AND 
HEAT TRANSFER 


boundary conditions must be applied is not known 
a priori. In order to make the equations converge 
satisfactorily, a rather compounded iteration pro- 
cedure is used. The transformed conservation 
equations and the expression for the heat transfer 
will be developed in the present section. The 
method of iteration and the numerical solution 
will be discussed in t he appendix. 

We first define a similarity variable for Ihe case 
of a sphere as 


v= 


V 


2 1/2 (uJL) 

PsPx 



( 40 ) 


Wo then express the streamwise velocity gradient 
•Mi and the normal velocity component r as follows: 


M| = -\ 


o 



(41) 


(42) 

At the stagnation point of the ax isy mine trie body, 
the continuity equation (18) is automatically 
satisfied b}^ the particular definition of J(y). The 
momentum, energy, and diffusion equations, equa- 
tions (19), (20), and (21), arc transformed as 
follows with the assumption of a constant Prandtl 
Dimmer: 


Mon entmn 


We shall, in the present section, solve the gov- 
erning equations, (18) through (21), with their 
boundary conditions for the velocity, enthalpy, 
temperature, atom concentration distributions, 
and shock-layer thickness. As a consequence of 
the solution, we shall then obtain the heat transfer 
to the stagnation point of spheres. 

The principles involved in solving the equations 
are rather straightforward. A boundary-layer- 
type affine transformation will be performed on 
the governing equations (18) through (21) and on 
their boundary conditions. The simultaneous 
ordinary differential equations thus obtained will 
then be solved numerically by the use of a digital 
computer. Though the principles involved are 
rather simple, the actual solution of the equations 
is quite complicated. The equations are strongly 
nonlinear and coupled. Moreover, the shock- 
wavi 1 detachment distance at which one set of the 


Enei gy 
Diffusion 




= 0 



(48) 


(44) 


(45) 


The reaction rate W/p is given by equation (39). 
The equation of state (22) and the expression of 
cntlt dpy given by equation (33) constitute two 
addi ional equations. The boundary conditions 
(23) through (28) become: at rj ~ 0 

/=/' = () ( 40 ) 
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at r) — ri 


f~ (50) 

\2 

//= 1 ( 51 ) 

C=i) (52) 

Thus there are five ('([nations, equations (43) 
through (45) and equations (22) and (3d), and five 
dependent variables /, p, T, (\ and //. W e nor- 
mallv need seven boundary eouditions for the 
present equations (43) through (45). There are, 
however, eight boundary conditions in the present 
analysis. The extra boundary condition is needed 
because t lie quantity r} s is an addition unknown. 
Now we see that the equations (43) through (52) 
and the equations (22) and (33) together constitute 
a consistent svstem lor a boundary value problem. 
The solution of the transformed equations will be 
obtained numerically with the use ot a digital 
computer. 

As was mentioned in the Introduction, the 
numerical integration of the equations requires 
considerable care because of the split boundary 
conditions and, most of all, the nonequilibrium 
chemical effect . In order not to complicate further 
the numerical work, here it is assumed that 

"U=l (53) 

for the entire stagnation-point shock layer. As- 
sumption (53) should not introduce any undue 
error into the beat-transfer results since it is the 
value of pg in the order of that near the shock 
which controls the heat transfer (see ref. 15). 
The degree of approximation involved in the 
present analysis due to assumption (53) will 
probably be of the same order as that involved in 
the hypersonic boundary-layer work of reference 
21 with the similar assumption. The analysis of 
reference 21 has been shown to be satisfactory 



for engineering purposes in comparison with the 
exact boundary-layer solution of reference 7. 

The properties of free-stream air used in the 
numerical analysis are obtained from reference 22. 
The IVandtl number is assumed to be 0.72 in the 
present- solution. 

The heat transfer is obtained from the solution 
of the equations as 

#SZ,— 2 I/4 \ PxMa 

RESULTS AND DISCUSSION 

Figures 2, 3, and 4 show three of the typical 
solutions of the conservation equations. All the 
■figures are for a flight speed of 20,000 ft /see, a 
nose radius of 1 foot, and a noncat alytie wall 
temperature of 1500 K. 

As was shown in reference 3, the vortieity inter- 
action which accompanies the thickening of the 
boundary layer begins to have a first-order effect- 
on the skin friction and heat transfer at the 
Knudsen number range corresponding to the 
present flight at 250,000 feet altitude. The 
velocity and the enthalpy profiles given in figure 
3 show the thickened boundary layer. Figure 4 
shows that the entire shock layer becomes quite 
viscous at a flight altitude of 280,000 feet. 1 he 
enthalpy profile shows that the effect of the cooled 
wall is felt by most of* the shock layer at this flight 
condition. 

The chemical state of the shock layer can be 
seen from the temperature and the atom mass- 
fraction profiles of figures 2, 3, and 4. In the 
conventional boundary-layer analysis, it has been 
assumed that the air is dissociated to its adiabatic 
equilibrium state at the edge of the boundary 
layer. The broken lines in figures 2, 3, and 4 
show the temperatures and the atom mass fractions 
which would be attained if the air were to reach 
equilibrium before the cooling eflect of the wall is 
felt. A study of the temperature and the atom 
mass fraction profiles in the light of the adiabatic 
equilibrium values of temperature and mass 
fraction shows that the entire shock layer is con- 
siderably out of equilibrium already at an altitude 
of 230,000 feet. As will be seen later in figure 5, 
the nonequilibrium chemical reaction in the shock 
layer at an altitude of about 230,000 feet causes 
a first-order deviation in the heat transfer from 
that calculated with the conventional assumption 
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(a) Velocity, enthalpy, and density profiles. 

Figure 4. — Typical solutions of conservation equations 
for u a . — 20,000 feet per second, L=] foot, T w — 1,500° 
K, and alt itude= 280,000 feet, Re^ 1,015. 



.2 .4 .6 .8 1.0 


y/A 

(b) Temperature and atom, mass-fraction profiles. 
Figure 4. — Concluded. 


of equilibrium boundary-layer-edge conditions. 
It should be remembered that the vorticity inter- 
action also begins to have a first-order effect on 
heat transfer at the similar flight condition (about 
250,000 ft). At the altitude of 280,000 feet (fig. 
4), it is seen that most of the shock layer is nearly 
frozen. 

Figure 5 shows the effect, of the chemical reaction 
on stagnation point heat transfer as a function of 
the flight altitude. The ordinate (tf/?i?) 8L 
obtained, in the present analysis, by the approxi- 
mate relationship, 

(55) 

The use of the above equation is justified here by 
the following argument. Our experience with 
the chemically reacting boundary layers and the 
comparison of the analyses of incompressible and 
compressible viscous shock layers found in 
references 1 and 3, respectively, show that the 
heat-transfer results are not greatly affected by the 
detailed variation of the fluid properties through 
the flow field. The main effect, therefore, of the 
chemical reaction on heat transfer to a noneatalytie 
wall appears through t he driving potential (h s — /<■„)« 

Figure 5 contains results for u m -~ 20,000 ft /sec 
and a result for TO, 000 ft /sec. Let us first 

consider the results for ?/„= 20,000 ft /sec. The 
solid line shows the heat-transfer ratio (q/^k) sl fur 



Figure 5. — Effect of chemical reaction on heat transfer to 
noneatalytie wall; u^ — 26,000 ft/sec unless otherwise 
specified. 
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the nose radius of 1 foot and the wall temperature 
of 1500° K. It is first seen from the solid line that 
the shock-layer dissociative chemical reaction 
reduces the heat transfer by the maximum relative 
amount when the flight altitude is about 210,000 
feet . At this altitude, the air dissociates essentially 
to its equilibrium state at the edge of the boundary 
layer. The boundary-layer reaction, which is a 
recombination process, on the other hand, is al- 
most completely frozen at this flight condition as 
was seen in reference 9. The shock layer disso- 
ciative relaxation and the near frozen boundary- 
layer recombination process together reduce the 
heat transfer by the maximum relative amount at 
the altitude of 210,000 feet. The heat-transfer 
ratio increases as the altitude deviates in either 
direction from the optimum height of 210,000 
feet. The heat-transfer ratio increases at the 
lower altitudes as the recombination process in- 
creases in the boundary layer. At the higher 
altitudes, on the other hand, the heat-transfer 
ratio (q /q E ) SL is increased as a smaller portion of 
the total energy becomes associated with the 
chemical energy of dissociation. This is because 
the degree of dissociation in the shock layer be- 
comes smaller as the dissociative relaxation time 
increases with altitude. The relaxation time 
increases as the pressure is decreased and also as 
the cooling effect of tin* wall is felt by a greater 
portion of the shock layer at the increased Knudsen 
number (see fig. 4). It can be said from tin 1 pre- 
ceding discussion that the dissociative relaxation 
in the shock layer is the controlling chemical re- 
action for heat transfer at the altitudes above 
about 210,000 teet, whereas it is the boundary- 
layer recombination which controls the heat trans- 
fer at the lower altitudes. The broken line in 
figure 5 is for the same conditions as the solid line 
except that the body nose radius is 5 feet. It is 
seen that the increased nose radius displaces the 
(qIQk)s L eurve to a higher altitude. The displace- 
ment is greater in the dissociation-controlled 
region than in the recombination-controlled region 
of lower altitudes. This means that the variation 
of the nose radius influences the shock-layer 
reaction more than it does the boundary-layer 
reaction. The main reason for this lies in the fact 
that the shock-layer thickness is approximately 
proportional to the nose radius, whereas the 
houndary-layer thickness is proportional to the 


square, root of the radius. It is seen that the air 
practically dissociates to the equilibrium state at 
the edge of the boundary layer at the approximate 
altitude of 240,000 feet when the nose radius is 
5 feet. 

The effects of the wall temperature on the heat- 
transfer ratio ( q!q E ) SL are shown in figure 5 by the 
filled symbols. At the altitude of 250,000 feet 
and for the nose radius of 1 foot, it is seen that the 
wall temperature has little effect on the dissocia- 
tive process in the shock layer. It is because the 
domain of influence of the wall temperature is 
still rather confined to the regions near the wall 
as seen by the enthalpy profile of figure 3. As 
the altitude and therefore the Knudsen number 
is increased, however, it is expected that the effect 
of the wall temperature will be felt by a greater 
portion of the shock layer (fig. 4) and, therefore, 
the effect of the wall temperature on the chemical 
reaction will become greater. The wall tempera- 
ture is seen to have a considerable effect on the 
heat t -ansfer through its effect on the boundary- 
layer recombination process at the lower altitudes. 
The boundary-layer reaction is influenced by the 
wall temperature substantially and this point was 
brought out in reference 9. 

A shock-layer solution for u w =30,000 ft/sec is 
shown in figure 5 for* a nose radius of 1 foot, wall 
temperature of 1500° K, and altitude of 250,000 
feet, ft is seen that the relative effect of the 
chemical miction on heat transfer is not too 
sensiti ve to the speed in the present (light regime. 

Before leaving figure 5, it may be important to 
note that the chemical reaction of dissociation and 
recombination in the stagnation-point shock layer 
ceases to be a problem above altitudes of about 
300,000 feet. 

We lave, so far, mainly discussed the effect of 
the chemical reaction on heat transfer. We have 
seen that the thickening of the boundary layer 
which accompanies the rarefaction of air has a 
large effect on the shock-layer chemical reaction. 
Beside* its effect on the chemical reaction, the 
thickei ing ol the boundary layer increases the 
heat ti ansfer directly. In figure 6, we shall see 
the simultaneous effect of the chemical reaction 
and the air rarefaction on the heat transfer. The 
curve of (qlq E ) SL versus altitude showing the 
chemic d reaction effects alone is transferred from 
figure 5 as the dashed line curve on figure 6. This 
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Figure 6. — Simultaneous effect of chemical reaction and 
air rarefaction on heat transfer to noncatalytic wall; 
w w = 26,000 ft/sec. 

ratio (<l!qu) SL shows only the chemical reaction 
effects because rarefaction effects occur in both 
the numerator and denominator. When the 
heating rate q SL is normalized by the beating rate, 
q HLyK} obtained from conventional equilibrium 
boundary-layer theory, t he ratio q &h !q B L^K will show 
the combined nonequilibriuni chemical and rare- 
faction effects since conventional equilibrium 
bounday-layer theory considers neither phenom- 
enon. This ratio is shown in figure 6 for two nose 
radii. The abscissa is again the flight altitude. 
From the present heat-transfer relationship (eq. 
(54)), and that given in reference 7, it is obtained 
that 

qsL 1 . 7 o () II J 

=— — - Y — fj — (oh) 

qiiL , E y /3 0 1 11 u\ E 

The streamwise velocity gradient at the edge of 
the boundary layer, (3 0) is obtained from the equa- 
tion given by reference 23, as 

0o=V«(2— e) 

where the density ratio e is here calculated by 
assuming that the air is in equilibrium immediate- 
ly behind the shock. 

At the lower altitudes, the classical boundary- 
layer theory is essentially valid and, therefore, 
the heat-transfer ratio qsL/qnL ,e is the same as the 
ratio (i qfqjs) SL * At the higher altitudes, figure 6 
shows that the heat-transfer ratio qsi.i Q bl,e be- 
comes higher than the ratio {q!qs) SL us the result 
of the direct effect of the air rarefaction (vorticity 


interaction) on heat transfer. On the solid curve 
corresponding to the 1-foot nose radius we note 
that at about 270,000 feet, the heat transfer is the 
same as that which would be calculated from con- 
ventional equilibrium boundary-layer theory. 
The reason is that the non-equilibrium chemical 
and air rarefaction effects which cause the heat 
transfer to deviate (in opposite directions) from 
the equilibrium boundary-layer value are just in 
balance and counteract one another at that point. 
At the altitudes below 270,000 feet, it is seen that 
the effect of the chemical reaction dominates over 
the direct effect of the air rarefaction, whereas the 
reverse is true at the higher altitudes. A similar 
trend occurs when the nose radius is 5 feet but at 
higher altitudes as shown in the figure. 

Finally, let us consider the simultaneous effect 
of the chemical reaction and the air rarefaction on 
the shock-layer thickness. The solid line and the 
filled points in figure 7 show the results of the 
present analysis. Figure 7 also contains the 
shock-layer thicknesses calculated by other in- 
vestigators for different conditions, which are 
shown here for the purpose of comparison. The 
shock-layer thicknesses based on the conventional 
inviscid flow theory are shown for € = 0.1 1 and for 
e based on the assumption of an equilibrium gas 
behind the shock. These values are obtained 
from reference 23. The shock-layer thickness cal- 
culated in reference 3 for the viscous compressible 
but chemically inert gas is also shown in t he figure. 
It is seen from the results of reference 3 that the 
shock layer becomes thinner as the air becomes 
rarefied and the effect of the cooled wall is felt by 
the greater portion of the shock layer. The result. 
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of the present analysis, on the other hand, shows 
that the shock-layer thickness increases as the 
air becomes rarefied. The explanation of this 
phenomenon is as follows. As the air is being 
rarefied, the shock-layer thickness is affected in 
t wo opposite ways. First, the cooling effect of the 
wall is felt by the greater portion of the shock 
layer and the shock layer tends to become thinner. 
With the rarefaction, however, the dissociative 
chemical reaction rate becomes slower and the 
degree of dissociation becomes smaller. The 
shock-layer thickness therefore tends to increase 
with the rarefaction of air. The over-all effect of 
the air rarefaction is that the chemical effect 
dominates over the cooled wall effect and the 
shock layer becomes thicker as the Reynolds num- 
ber is decreased in the present flight regime. 

CONCLUDING REMARKS 

The nonequilibrium chemical reaction of dis- 
sociation and recombination is studied theoreti- 
cally for air in the viscous shock layer at the 
stagnation region of axisymmetrie bodies. The 
flight regime considered is for speeds near satellite 
speed and for altitudes between 200,000 and 
300,000 feet. 

As altitude is increased from 200,000 to 300,000 
feet, the flow in the shock layer varies from one 
which supports a well-defined boundary layer to 
one characterized by a thick viscous layer. Tin* 
chemical reaction rate, at the same time, varies 
from near equilibrium to completely frozen. 

The nonequilibrium chemical reaction is found 
to produce a significant reduction in the heat 
transfer to noncat alytic walls at certain altitudes. 
For a wall temperature of 1500° K and nose 
radius of 1 foot the significant altitude lies 
between about 150,000 and 270,000 feet. The 
maximum reduction in heat transfer occurs at 
about 210,000 feet altitude and the heat transfer 
is then about one- third of the equilibrium value. 

1 1 is shown that the controlling chemical reaction 
for heat transfer at higher altitudes is the dis- 
sociative shock-layer reaction, whereas at lower 
altitudes it is the boundary-layer recombination 
process. 


Th k air at the boundary-layer edge is found to 
be ir a state of chemical equilibrium below 
altitudes between 210,000 and 240,000 feet for 
nose -adii from 1 to 5 feet. Most of t lie shock 
layer, however, is found to be substantially out of 
equilibrium throughout the flight regime con- 
sidered. It can be said from this that the radia- 
tion characteristics from the shock layer may he 
considerably different from that calculated with 
the assumption of an equilibrium shock layer. 
Also, the gas which surrounds a vehicle and the 
gas a the wake and trail of the vehicle probably 
contam substantially less amounts of dissociated 
radicals than they would if the stagnation shock 
layer were in equilibrium. 

Th«‘ relative effect of chemical reaction on heat 
transfer is found to be rather insensitive to flight 
speed. 

Nose radius is seen to affect the dissociative 
chemical reaction in the shock layer more than 
it does the recombination process in the boundary 
layer. The wall temperature, on the other hand, 
influences the boundary -layer reaction at lower 
alt it u les most significantly, and the shock-layer 
dissociative reaction at extremely high altitudes 
mode utely. At intermediate altitudes, however, 
the effect of wall temperature on chemical ruction 
is negligible. 

As the altitude is increased, rarefaction of the 
air is .'ound to produce two opposite effects on the 
shock layer thickness. First, the effect of the 
cooler wall is felt by a greater portion of the shock 
layer md the shock layer tends to become thinner. 
The j arefaction, at the same time, freezes the 
dissociative chemical reaction which tends to 
increase the thickness. The overall effect of the 
air n refaction on the shoek-layer thickness is 
that t ic chemical effect dominates over the cooled 
wall effect and the shock layer becomes thicker as 
the ai * is rarefied. Finally, the chemical reaction 
kineti *s at the temperatures above about 7000° K 
are n >t well known at the present time. The 
presei t results, therefore, should he considered 
with 1 his fact in mind. 

Ames Teseakch Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., May 8, 1!)61 



APPENDIX A 


ITERATION PROCEDURE AND NUMERICAL SOLUTION OF CONSERVATION EQUATIONS 


Some of the basic principles programmed into 
the digital computer to solve equations (43), (44), 
and (45) arc discussed in this appendix. 

A conventional, straightforward iteration tech- 
nique is used to integrate the momentum and 
the energy equations. The solution of the 
momentum equation begins with the guessing of 
a function p( 77). The two known boundary 
conditions (eq. (46)) are then applied, and a third 
boundary condition /" (0) is guessed at the wall. 
The forward integration then begins and con- 
tinues until one of the two boundary conditions, 
say equation (49), is satisfied. At this point, the 
other boundary condition (50) is checked to see 
if it is also satisfied. If the boundary condition 
(50) is not satisfied, then another value of /" (0) 
is guessed and the equation is again integrated. 
The iteration thus continues on the boundary 
condition (50), with the aid of /" (0), until all the 
boundary conditions are satisfied. The solution 
of the momentum equation gives Vt as well as 
f(v) and its derivatives. The method of solution 
of the energy equation is quite similar. For a 
given wall temperature, C w is guessed and this 
yields the boundary condition II W with the aid 
of equation (47). Another boundary condition 
II' (0) is guessed at the wall and the forward 
integration begins. The integration continues 
to Vs at which point a check is made to see if the 
boundary condition (51) is satisfied. If not, a 
new guess is made on TV (0) and the integration 
is repeated. The iteration is continued until 
all the boundary conditions are satisfied. 

The diffusion equation (45) is not amenable to 
such straightforward iteration procedure. There 
are essentially two categories of difficulties one 
must face in order to solve the equation suc- 
cessfully. 

The first category of difficulty arises with the 
increasing magnitude of the reaction rate term 
(L/u m ) (W/p). As the magnitude of the chemical 
reaction term (L/u m ) ( W/p ) is increased, the 
equation fails to converge unless almost exactly 


the correct value of C w lias been guessed and 
H(v) has been obtained. For most of the reason- 
able values guessed, not only does the equation 
fail to converge but it diverges even before the 
integration can reach v s . To circumvent this 
difficulty, the following procedure was used in 
the present analysis. The equation to be solved is 


_ 1 _ 

Pr 


< 7 "+ 2 fC' = - 


V2 



(45) 


with the boundary conditions 
6 y , ( 0)=0 
CM = 0 


(48) 

(52) 


The main trouble here is caused by the nonlinear 
effect of the term (L/u^) (W/p). The energy and 
the diffusion equations are iterated together in the 
present endeavor. With/(iy) known from the so- 
lution of equation (43), diffusion equation (45) is 
linear for each iterative step if (W/p) is calculated 
from the values of T and C obtained from the pre- 
ceding iteration of the energy and diffusion equa- 
tions. Therefore, for the nth iteration step and 
for a given f(v) we shall linearize the diffusion 
equation as 


1 

Pr 


(\"+2fC n ' = - 


1 L 
v 2 u m 



r 


(Al) 


with tlie boundary conditions 


C n ' (())=() c n (v ,)= 0 


(A2) 


The symbol T represents a positive number be- 
tween 0 and 1, and is put in here to further insure 
the convergence. During the iteration for a par- 
ticular n, if the diffusion equation diverges before 
v reaches v* y the computer is programmed to reduce 
the effect of W/p by reducing T until the integra- 
tion can be accomplished satisfactorily. The re- 
sults for the reduced V are then used to calculate 
a new W/p, and the diffusion equation is rein- 

15 
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tegrated for a new F wliich is greater than the 
prereding one by a certain suitable amount. The 
solution of equation (Al) for the particular n is 
built up, in this manner, until F becomes unity. 
When F becomes unity, the simultaneous iteration 
of the energy and the diffusion equations resumes. 
The iteration of the two equations continues until 
they converge. The function p(rj) is then recal- 
culated and the momentum equation is solved with 
the new p(t?). The final solution of the problem 
is obtained when all three equations converge 
simultaneously. 

The second category of difficulty arises during 
the course of the iterations explained above. The 
difficulty is associated with the additional iteration 
required to satisfy the split boundary conditions 
(A2) of the diffusion equation. To circumvent 
this difficulty, we express the solution of equation 
(Al) by 

- B n C c , n (v) +C P , n (rj) (A3) 

The above expression is permissible because 
equation (Al) is linear. The symbol C e , n (y) 
represents the complementary function, whereas 
(’p.niv) represents the particular solution. The 
function C CtH , therefore, is a solution of the 
homogeneous equation 

A C f „"+2/C c ./=0 (A4) 


whereas (' Pt „ is a solution of the nonhomogeneous 
equation 


1 

Pr 


(\ 


'+2/fV 



(AS) 


We obtain solutions C Ct „ and C p , n to satisfy the 
follov/ing boundary conditions, respectively, all 
at 77 - =0, 


^ v „ (o) = 1 9 

(A 6) 

<\,n («) = (> J 

<V.(o)=o1 

(A7) 

(\.n (0) = 0 J 

jj ( p,n (^x) 

n ~ t'ZM 

(AS) 


Now it is seen that (A3) is the complete solution 
of equation (Al) which satisfies the boundary 
conditions (A2). Let us now examine equations 
(A4) and (A5). An examination of equation (A4) 
in the light of the boundary conditions (AG) shows 
that the unique solution of the homogeneous 
equation is C Ctn {y) = l. The complete solution for 
the nth iterated equation (Al) can therefore he 
written as 


( y n {y) = — C Ptn (rf t ) + C Pt n (rj) (A9) 

The integration which must be performed for the 
diffusion equation with a digital computer is now 
reduced to the integration of equation (A5) with 
the boundary conditions (A7). The boundary 
conditions (A7) are both given at t? = 0 and, there- 
fore, K) iteration is involved in the solution of 
equation (AS). 
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